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ABSTRACT 

The dilation of a Euclidean graph is defined as the ratio of 
distance in the graph divided by distance in R''. In this pa- 
per we consider the problem of positioning the root of a star 
such that the dilation of the resulting star is minimal. We 
present a deterministic 0(n log n)-time algorithm for evalu- 
ating the dilation of a given star; a randomized 0{n\ogn) 
expected-time algorithm for finding an optimal center in 
R'*; and for the case d = 2, a randomized 0(n2"'"' log^ n) 
expected-time algorithm for finding an optimal center among 
the input points. 

Categories and Subject Descriptors 

F.2.2 [Theory of Computation]: Analysis of Algorithms 
and Problem Complexity — Nonnumerical Algorithms and Prob- 
lems; G.1.6 [Mathematics of Computing]: Numerical 
Analysis — Optimization 

General Terms 

Algorithms, performance, theory 

Keywords 

Facility location, quasiconvex programming, dilation, stretch 
factor, Euclidean graphs 



-1 1. INTRODUCTION 



A star is a graph with exactly one internal vertex, called 
its center; it has edges from the center to every external 
vertex, and no other edges. The dilation between any pair 
of vertices a and 6 in a Euclidean graph is defined [5] as 
the cost of the shortest path from a to b, divided by the 
Euclidean distance \ab\ . The dilation of a graph is defined as 
the maximum dilation over all pairs of vertices. We consider 
the following problem: for any fixed d > 2, given a set 
1/ C R'* of n points, construct a star with center c G R'* 
and leaves V such that the star has minimal dilation. We 
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present algorithms for the case when c is constrained to be 
one of the input points, as well as the unconstrained case 
when c may be any point in R'*. 

We present algorithms to evaluate the dilation of a given 
star in 0{n\ogn) time; find the optimal center c £ R'* for 
a given set of external vertices in 0(n log n) expected time; 
and for a set of vertices 1/ C R'^, select an optimal center 
c G V in expected 0(n 2°'''"-' log'^ n) time. The evaluation al- 
gorithm's task is to find the single pair of points that define 
the graph's dilation. It works by identifying 0(n) pairs of 
points with the potential to be that pair in O(nlogn) time, 
and evaluating the dilation for each of those 0{n) point 
pairs. Our algorithm for the unconstrained optimization al- 
gorithm formulates the problem as a quasiconvex program, 
and uses the evaluation algorithm as a component of Timo- 
thy Chan's randomized optimization framework [3] to arrive 
at the solution in 0(n log n) expected time. The algorithm 
for the constrained case works by repeatedly selecting a ran- 
dom vertex, evaluating the dilation that would result from 
using that vertex as the center, computing the region R of 
space that could contain a center yielding an even lower di- 
lation, and discarding all the vertices outside R. This proce- 
dure is iterated an expected O(logn) times, and computing 
a description of R takes 0(n 2"'-"^ log n) time, resulting in 
an 0{n 2°'''"-' log'^ n) expected-time algorithm. 

Finding a minimum dilation star can be viewed as an in- 
stance of the classic facility location problem [4] : given a set 
of demand sites and supply sites, find a center minimizing 
a particular objective function. In our case the objective 
function is defined as the maximum dilation between any 
pair of points. This formulation could e.g. represent the 
problem of deciding where to position an airline hub, such 
that flights with layovers are not unreasonably longer than 
direct flights. 

Narasimhan and Smid [8] considered the related problem 
of computing dilation approximately. They present algo- 
rithms which can compute an e-approximation of the dila- 
tion of a path in O(nlogn) time, a cycle in 0(n log n) time, 
and a tree in 0(n log'^ n) time, for any flxed d > 1. Agarwal 
et al. [1] considered the problem of computing the dilation 
of a polygonal curve. They present algorithms to compute 
the exact dilation of a polygonal curve in R'^ in 0{nlogn) 
expected time, or in R"^ in subquadratic time. Langerman 
et al. [7] presented algorithms which compute the dilation of 
a planar polygonal curve in expected 0(n log n) time, a pla- 
nar tree in 0{n log^ n) time, or a planar cycle in 0(n 2 log n) 
time. The evaluation algorithm presented in this paper com- 



putes the exact dilation of a star in 0(n log n) time for any 
fixed d > 2. 

2. THE EVALUATION PROBLEM 

Formally, the evaluation problem is as follows: given a set 
V oin points in R** and a point c £ R'', compute the dilation 
of the Euclidean star with center c and leaves V . 

For any pair of vertices a and b, a star contains exactly 
one path from a to b: the path {a,c,b). So the dilation 
between a and b is 



Ac{a,b) 



\ac\ + \ cb\ 
\ab\ 



where \ac\ designates the Euclidean distance from a to c. 
The dilation of the star is 

Ac — max Ac(a,b) . 

a,beV 

Define a to be one of the points in V for which Ac{a,a) 
is greatest (there may be more than one pair of points with 
dilation Ac). As noted in the introduction, one approach 
to solving this problem is to reason that for any a £ V, a 
could be any of the n — 1 points in 1/ \ {a}. Then Ac may 
be computed by an algorithm that evaluates Ac (a, b) for all 
O(n^) pairs of points. As we will see, this approach is overly 
conservative. In this section we present an algorithm that 
identifies 0(n) candidate pairs of points, such that the pair 
(a, a) for which Ac = Ac(a, a) is guaranteed to be one of the 
candidates. The algorithm takes 0(n log n) time to generate 
the candidate list, and computing the largest dilation among 
0(n) pairs of points takes 0{n) time. The point pairs are 
identified using two techniques, each generating 0(n) pairs. 
One technique generates a list of 0(n) pairs of points guar- 
anteed to contain (o, a) when Ac is high; the other generates 
0(n) pairs guaranteed to contain (a, a) when Ac is low. 

2.1 The high dilation case 

One of our heuristics for finding a is to identify a's k 
nearest neighbors, for a constant k £ 0{1). This heuristic 
is guaranteed to identify a when Ac — Ac (a, a) is high; 
formally, we will prove that if Ac > F for a constant F > 3, 
then a is one of a's k nearest neighbors, for a constant k. 
The k nearest neighbors of every a £ V may be reported in 
0(fcn log n) time using the algorithm of Vaidya [10], so the 
process of identifying these 0(n) point pairs takes 0{n log n) 
time. 

Lemma 1. Suppose a £ V and Ac = Ac(a,a) > F for a 
constant F > 3. Then \ac\ > (j!)|aa| for a constant (ji > 1 that 
depends only on F. 

Proof. By definition. 



F < Ac = ' < + + < 1 + 2M 



ac + \ca\ \ac\ 



I I r- 1, . 

ac > \aa\ 

I I — oil 



Therefore if F > 3, = > 1. 



□ 



Lemma 2. Suppose that Ac — Ac{a,a) > F for a con- 
stant F > 3, let K be the d-sphere centered on a with ra- 
dius \aa\, and u,v £V he two input points inside k. Then 




Figure 1: k partitioned into small d-cubes for d = 2. 

there is a constant 7 > that depends only on F, such that 
\uv\ > 7|aa|. 

Proof. By Lemma0 \ac\ — \aa\ > (1 — l/<^)|ac|, so 



3|ac| lacl + \ac\ . . / s 

\aa\ \aa\ 

^ \uc\ + \cv\ ^ 2(|ac|-|aa|) 
~ \uv\ ~ \uv\ 

^ 2(l-l/,^)|ac| 



It follows that > |(1 — l/<f>)\aa\. So \uv\ > 7|aa| for 



7 = f (1- > 1 so 7 > 0. 



□ 



Lemma 3. If Ac > F and Ac — Ac{a,a), then a ts one 
of a's k nearest neighbors, for a constant k depending only 
on F and the dimension. 

Proof. As shown in Lemma|21 the distance between any 
u and V is proportional to the radius = |aa| of k. Let 
(T = ^/y/d and partition n into small d-cubes with sides of 
length (Tr„ (see Figure 0. Under such an arrangement no 
cube contains more than one input point, k, can contain no 
more than [(§7^)''] ~ [(J-)'*! such cubes, which is constant 
in n. So we may conclude that if Ac = Ac(a, a) then a is 
one of a's k nearest neighbors, for a constant k = [(-)'']. 

□ 

2.2 The low dilation case 

Our other heuristic for finding a is to identify points whose 
distance from c is similar to that of a. This may be done 
by sorting V by distance from c, then for each a, identifying 
the points whose indices in the sorted sequence are within 
/ of a's, for a constant I. We will show that this heuristic 
is guaranteed to identify a when Ac = Ac(a,d) is low, i.e., 
that if Ac < F for a constant F, and Vs is the set V sorted 
by distance from c, then for any Vi € Vs, if Vj = Vi and 
Ac = Ac{yi,Vj) then \i — j\ < I, for a constant I. We 
prove this by showing that space may be partitioned into 
d-dimensional annuli with exponentially growing radii, all 
centered on c, such that each annulus contains 0(1) points 
from V, and for any a £ V, a lies within one of 0(1) adjacent 
annuli. We define an annulus Ac,i in terms of its inner and 
outer radii, and p'"^^, respectively. 



A.,, = {p G R'' 1 < \pc\ < p'+'} 
where the constant p is defined as 



A, + 1 



A, 



Lemma 4. If a £ Ac,i and a G Ac,i+k, then k <2. 

Proof. Suppose for the sake of contradiction that fc > 3. 
We now establish an upper bound on Ac(a, a) by positioning 
a and a such that dilation is maximized. This happens when 
a and a are colinear and as close as possible. We may assume 
without loss of generality \ac\ < \ac\; so \ac\ = p'"*"^ and 
\ac\ = p*"*"^, hence 



Ada, 



(A, 



- Ae+ 1 



Ac + 1 _ 1 



Ac + 1- 



= Ac 



which is a contradiction; so the assumption fc > 3 must be 
false. 

□ 

Lemma 5. There exists a function 6min such that if a 
and b lie m the same annulus, Zacb > SM/iv(Ac), and a 
value Ao such that for all Ac < Ao, Sm/jv(Ac) > n/2. 

Proof. We wish to establish a lower bound on Zacb, and 
so we consider the choice of a and b that makes Zacb as acute 
as possible. For a fixed Ac (a, 6), we can do this by making 
||ac| — |c6|| as large as possible. Without loss of generality, 
assume \cb\ = p\ac\; then 



Ac > 



\ac\ + p\ac\ 
\ab\ 



Let 9 — Zacb; then by the law of cosines, 

\ac\ + p\ac\ 



Ac > 



\/|acp + p-^lacp — 2p|acp cos^ 
1 + P 



2p cos 6 



Solving for 



so we let 9 a 



2p 2pAi 



KAc) = cos-i(f + i 



2p 



follows from the convergence of the limit 
lim 9min(Ac) = vr . 



(1+pr 

2pA^ 



The result 



□ 



Lemma 6. Each Ac,i contains 0{1) points from V . 

Proof. Lemma |K| shows that if a and b axe in the same 
annulus and 6 = Zacb, then 



Thus each annulus contains 0{9q points; do and d are 
constant with respect to n, so this quantity is as well. 

□ 

Lemma 7. // Ac < F, Vs is the set V sorted by distance 
from c, Vs and vt G Vs, and Ac — Ac{vs,vt), then \s — t\ <l 
for a constant I depending only on F and d. 

Proof. Lemma |1| shows that if a G Ac,i and a G Acj, 
then \i — j\ < 2. The points in contiguous annuli will be 
contiguous in Vs- So a must be within I ranks of a in T^s, 
where I is the number of points that may lie in two annuli. 
As shown in Lemma El I is constant with respect to n. 

□ 

2.3 The evaluation algorithm 

Theorem 1. Given a set V C of n points and center 
c G R"^, it is possible to compute the dilation Ac for the star 
with center c and leaves V in 0(n log n) time, provided d is 
constant. 

Proof. We generate a list of 0(n) pairs of points as fol- 
lows. Compute the k nearest neighbors for the set V in 
0(n log n) time, for the constant k described in Subsec- 
tion 12.11 Then for each a £ V, append the point pairs 
consisting of a paired with each of a's k nearest neighbors 
to the list. Next, create Vs by sorting the points in V by 
distance from c in 0(n log n) time. For each Vi G Vs, add 
the pair {vi,Vj) to the list, for every j such that \i — j\ < I, 
for the constant I described in Subsection l2.2l The resulting 
list has (kn + In) G 0(n) elements. As shown in Lemmas El 
and[7| the pair defining the dilation of the star, (a, a) such 
that Ac = Ac(a,a), is certain to be present in the list. So 
we can compute Ac by evaluating the dilation of each of the 
0{n) pairs of points, and returning the maximum. 

□ 

3. THE UNCONSTRAINED 
OPTIMIZATION PROBLEM 

We now turn to the following problem: given a set of n 
points V C R'', find a point c, not necessarily belonging to 
V, such that the star with center c and leaves V has minimal 
dilation. 

3.1 Reduction to a quasiconvex program on 
O(n^) contraints 

A function / : R*^ i-^ R is quasiconvex when its level 
sets /-^ = {a; G R"^ I f{x) < A} are convex. Quasicon- 
vex programming [6] is a generalization of linear program- 
ming: a quasiconvex program is a finite set S of quasiconvex 
functions (constraints); the solution to a quasiconvex pro- 
gram is the input a; G R"* minimizing the objective function 
fs = max/^gg fi{x). Amenta et al. showed [2] that if a so- 
lution to a constant-sized subset of S may be solved in 0(1) 
time, then the program may be solved in 0{\S\) expected 
time. 

Lemma 8. The star center c minimizing Ac may be com- 
puted by solving a quasiconvex program with 0(n^) con- 
straints. 

rid , 



Proof. We introduce one contraint /^j 
each pair of points Vi,Vj G V: 



for 



\vix\ -\- IxVil 




Figure 2: An example level set /^^ . 

Hence the solution x is precisely the point c for which the 
dilation is minimized for all pairs of input points. 

The denominator of fi^j is the distance between Vi and 
Vj, which is constant with respect to x. So fi,j{x) is pro- 
portional to the sum of the distances of x from Vi and Vj , or 
equivalently each fi,j defines an ellipsoid with foci Vi and Vj. 
So any level set ff^ is elliptical (see Figure |5J, and hence 
convex; so each fij is quasiconvex. 

□ 

3.2 The optimization algorithm 

In this section we show how to solve the unconstrained op- 
timization problem in O(nlogn) expected time. We adapt 
Timothy Chan's randomized optimization technique [3] to 
our problem. Chan states his result in terms of LP-type 
problems, but following [6], we rephrase it in terms of qua- 
siconvex programming; 

Lemma 9 (Chan). Let Q be a space of quasiconvex func- 
tions, V be a space of input values, and / : 2^^ 2^ map 
sets of input values to sets of functions in Q. Further, sup- 
pose that V , f , and S satisfy the following properties: 

• There exists a constant-time subroutine for solving qua- 
siconvex programs of the form f{B) for any B C V 
with \B\ G 0(1). 

• There exists a decision algorithm that takes as input a 
set P C V and a pair (A, x), and returns yes if and only 
if q{x) < X for all q G f{P)- The running time of the 
decision algorithm is bounded by D{\P\), where there 
exists a constant e > such that D{n)/n'^ is monotone 
increasing. 

• There are constants s and r such that, for any input set 
P G V , we can find in time at most D{\P\) a collection 
of sets Pi, < i < r, each of size at most s\P\, for 
which /(P) =U,/(PO- 

Then for any P C P we can solve the quasiconvex program 
f{P), where \P\ — n, in randomized expected time 0{D{n)). 

We now show how to apply this result to the problem at 
hand. 

Theorem 2. Given a set V of n points, it is possible to 
compute a point c that admits the minimal- dilation star with 
center c and leaves V in expected 0{n log n) time. 

Proof. Lemma |S| shows how to formulate this problem 
as a LP-type program with quasiconvex constraints. Sec- 
tion |5| provides an 0{n log n) time algorithm to evaluate the 



dilation of a star with n leaves; so our decision algorithm 
can invoke that algorithm, then compare the result to A, to 
decide whether q{x) < A for any A. The resulting decision 
algorithm has running time D(n) G 0(n log n). 

A set of input points P may be divided into r subsets, each 
with size [s|P|], as follows. To ensure correctness, every 
pairing of points must be mutually present in at least one 
Pi. So we partition P into three arbitrary disjoint subsets 
Qi,Q2, and Qa, each of size \Qi\ > L|P|/3J, and form three 
subsets Pi = P \ gi, P2 = P \ Q2, and Pa = P \ Qs- Then 
any pair of points present in P is also mutually present in 
some Pi. Each Pi has size \Pi\ < ||P|, so our approach 
fits the requirements of Chan's framework with parameters 
r = 3 and s = | . Consequently it may be solved in expected 
OiD{n)) = 0(n log n) time. 

□ 

4. THE CONSTRAINED OPTIMIZATION 
PROBLEM 

Finally we consider the following problem: given a set 
y C of n points, find a point c € V such that the star 
with center c and leaves V \ {c} has minimal dilation. 

At a high level, our algorithm is as follows; 

CONSTRAINED-DILATIONCl^); 

1. Let C ^ y 

2. Repeat 

(a) Let c be a random point from C. 

(b) Compute Ac, the dilation of the star with center c 
and leaves V \ {c}, using the algorithm presented 
in Section 121 

(c) Compute the region R which may contain a center 
c' that admits a dilation A^.' < Ac. 

(d) C^CnR. 

(e) If C = 0, return c. 

In each iteration of the loop, any a £ C will be removed 
with probability 1/2; so the expected number of iterations 
is log n. The remainder of this section will describe a rep- 
resentation of R that may be constructed and queried ef- 
ficiently. We will present a simple data structure that can 
be constructed in 0(n 2"'"' log n) time and answer member- 
ship queries in O(logn) time, which makes the overall ex- 
pected running time of CONSTRAINED - DILATION 
0(n2"("^ log^n). 

The region R is the intersection of a number of ellipses. 
We will use an argument similar to that of Section |5| to 
show that R may be accurately represented by 0{n) ellipses, 
which may be identified in 0(n log n) time. Our algorithm 
needs to be able to decide whether given points lie inside, 
or outside, the intersection of those 0{n) ellipses. We will 
show that this intersection can be described as a sequence of 
arcs. A straightforward divide-and-conquer algo- 
rithm can construct this sequence in 0(71 2"'-"^ logn) time, 
and binary searches in the sequence can answer membership 
queries in O(logn) time. 

Lemma 10. For any quantity Ac and V C R'^, the region 
R = {r \ Ar < Ac} ts defined by the intersection of ellipses. 




Figure 3: Arcs defining an intersection of ellipses. 

Proof. The proof of Lemma|H|contains a description for 
how to form a constraint fi,j for each pair of points in V. As 
discussed in that proof, the level set f^^ defines the region of 
JR"* for which Ac{vi,Vj) < X. So the region R is the interior 
of the intersection of these ellipses; formally, 

i? = int n f^f-^ 

(see Figure 121. 

□ 

Lemma 11. If Ac > F for a constant F > 3, then the 
region R is the intersection of 0{n) ellipses, which may be 
identified in 0(n log n) time. 

Proof. Let a and b be the foci of an ellipse contributing 
to the boundary of R and containing the boundary point c'; 
then Ac = A^' = Ac'(a,&). By Lemma |21 b must be one 
of the k nearest neighbors of a. Therefore, the only ellipses 
that contribute to the boundary of R are those found by 
reporting the k nearest neighbors of each input point. 

□ 

Lemma 12. //Ac < F for a constant F, then the region R 
is the intersection of 0{n) ellipses, which may be identified 
in 0(n log n) expected time. 

Proof. As in Lemma [TTI we will use the algorithm pre- 
sented in Section |5|to identify the relevant ellipses. However 
the algorithm presented in Subsection 12.21 expects a center c 
as input, and we are now computing c as output. We over- 
come this obstacle by computing cqpt, the solution to the 
unconstrained optimization problem discussed in Section |3 
We then use the approach of Subsection l2.2l simulating the 
annuli centered on c with the annuli centered on cqpt- We 
will argue that for any a £V we can use techniques similar 
to those in Subsection 12 . 21 to identify 0(1) points, such that 
a is one of those points. Let x = |ccopt|; then a may fall 
into one of three categories. It may be roughly x away from 
Copt, or closer, or farther. 

Let a £ V, and suppose p~^x < \acoPT\ < px. Then 
a can only lie in one of two annuli centered on copt, i.e. 
AcoPT,i or ^copT,i+ij for ^ — [logp ; so for any a we may 
pessimistically include all the 0(1) points in those annuli in 
the list of a candidates. 



Suppose \acopT\ > x; then shifting the annuli's center 
from c to COPT will cause a to cross at most one annulus' 
border. So if a e Ac,i then a G AcQp^.i-i U Ac^pT^i U 
AcoPT,i+i-- So for these points we include the points within 
21 ranks of a in Vs, rather than I ranks as in Subsection l2.2l 

Finally we consider points a such that |acoPT| < x. By 
the assumption that Ac^p^ < F there can only be 0(1) 
such points; so we can pessimistically include all 0(1) such 
points. 

So for each a £ V, we form the ellipses corresponding to 
a paired with the points in the two annuli whose radii are 
close to x; the points appearing within 21 ranks of a in the 
sorted sequence Vs; and the 0(1) points within x of copt- 

□ 

Definition 1. An arc ring is a pair {r,S) describing the 
convex region inside the intersection of ellipses, r is an ar- 
bitrary point inside the region. S = ((ei, 6'i), (e2, ^2), . • .) 
is a sequence of arcs, with each element {ei,6i) describing 
a range of angles 9i-i < 9 < 9i about r, for which the el- 
lipse Ci defines the boundary of the region (where 9o = 6^g^). 
We reguire that the angle boundaries 9i appear in ascending 
order, 9i — 0, and 6^s\ = Stt. 

Lemma 13. There exists an arc ring describing any R us- 
ing 0(n2"(")) arcs. 

Proof. Lemmas llll and lT^ show that for any value of Ac, 
R may be described as the intersection of 0{n) ellipses. As 
shown in [9], a description of the intersection of 0{n) ellipses 
forms a Davenport-Schinzel sequence of order four; such a 
sequence has length 0(n2"("'). Hence \S\ G 0(n2"(")). 

□ 

Lemma 14. An arc ring E describing R may be generated 
in 0(n 2°''"-' log n) time and queried in O(logn) time. 

Proof. E may be formed in 0{n 2°'"' log n) time using a 
simple divide-and-conquer algorithm resembling merge sort. 
First partition the set of ellipses into two subsets of equal 
size, and compute the intersections Si and S2 of the sub- 
sequences recursively. Then merge Si and S2, as in merge 
sort, by comparing the least elements of each subsequence, 
si and 82 respectively, computing an element to append to 
S = SiU S2, and modifying or deleting si and S2. This may 
be done by partitioning the element with a greater angular 
interval into two elements, so that si and S2 cover exactly 
the same interval. If si and S2 do not intersect over that 
interval, then append whichever element is closer to r to S; 
otherwise, append two elements corresponding to the closest 
ellipse on either side of the intersection point. The running 
time of this algorithm is given by the recurrence 

r(|S|) = 2T(lSl/2) + 0(|Sl) , 

which solves to T{n) G OdSI log \S\); so the running time is 
0(n2''("hog(n2"("')) = 0(n 2"'") log n). Once has been 
constructed, it is possible to determine whether a point a 
lies in R by computing the angle between r and a, then 
performing a binary search in S to find the pertinent arc, 
and computing whether a lies inside or outside that arc. 
This operation takes O(logn) time. 

□ 

Theorem 3. The optimal center c £ V may be found in 
0{n 2"'"' log^ n) expected time. 



Proof. The expected number of iterations of the outer 
loop of CONSTRAINED - DILATION is log n. The 
running time of each iteration is dominated by the con- 
struction of R, which can be achieved in 0(n 2°'"' log n) 
time. Hence the algorithm runs in 0(n 2«(">log2 n) ex- 
pected time. 

□ 
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